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Abstract 

Cubic invariants for two-dimensional Hamiltonian systems are investigated using the Jacobi 
geometrization procedure. This approach allows for a unified treatment of invariants at both fixed 
and arbitrary energy. In the geometric picture the invariant generally corresponds to a third rank 
Killing tensor, whose existence at a fixed energy value forces the metric to satisfy a nonlinear 
integrability condition expressed in terms of a Kahler potential. Further conditions, leading to 
a system of equations which is overdetermined except for singular cases, are added when the 
energy is arbitrary. As solutions to these equations we obtain several new superintegrable cases 
in addition to the previously known cases. We also discover a superintegrable case where the 
cubic invariant is of a new type which can be represented by an energy dependent linear invariant. 
A complete list of all known systems which admit a cubic invariant at arbitrary energy is given. 
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I Introduction 



The quest for integrable systems has a long history going back more than a century. However, to this day, 
no general systematic method exists for finding the invariants of a given system. An effective criterion to 
test if a linear invariant exists in a 2-dimensional system was found only recently in jlj| (compare also Q). 
In the absence of such general methods one can at least study and classify sufficiently simple integrable 
systems with certain simple types of invariants. In the present work, we focus on cubic invariants admitted 
by natural Hamiltonian systems in two-dimensional Euclidean space. 

Cubic invariants in two dimensions have been studied before by several workers. In 1935 Drach || carried 
out the first systematic study of such systems. However, using null variables without restricting them to be 
complex conjugate, Drach's Hamiltonians were not constrained to be the sum of the Euclidean space kinetic 
energy and a real potential. This is probably why his results have not often been compared to those of more 
recent studies such as that of Holt in 1982 ||]. The most complete classification of 2-dimensional systems 
admitting a cubic invariant, but not including the Drach systems, was later given by Hietarinta in 1987 [^). 
Most workers have dealt with arbitraray energy (strongly conserved) invariants, but fixed energy (weakly 
conserved) invariants have also been investigated to some extent (c.f. Q). In this paper we use the Jacobi 
geometrization method to represent the dynamics. That approach was used in |6| to give a unified treatment 
of quadratic invariants at fixed and arbitrary energy. We will show that such a unification is possible also 
for systems with cubic second invariants. 

The Jacobi formulation has the advantage of encoding the entire dynamics in a single geometric object, 
the Jacobi metric. The system orbits are simply the geodesies of the Jacobi geometry. To analyze the system 
one can utilize all the mathematical and computer algebra tools developed for Riemannian geometry. Other 
geometric methods which have been used for integrable systems can be found, e.g., in Q|. In particular we 
mention the projection method of Olshanetsky and Perelomov That formulation, however, is fundamen- 
tally different from the Jacobi geometrization. For example, the projection geometry itself is not sufficient to 
represent the dynamics. On the other hand the projection geometry is designed to be very simple, typically 
a space of constant curvature. A drawback is that there is no general algorithm to find a suitable projection. 

A polynomial invariant of the geodesic equations corresponds to a single geometric object on the con- 
figuration space, known as a Killing tensor. The vanishing of the Poisson commutator of the invariant and 
the Hamiltonian gives a set of conditions which are referred to as the Killing tensor equations. Because 
of Noether's theorem, such invariants correspond to certain dynamical symmetries of the geodesic equa- 
tions |^| . When dealing with cubic invariants, in contrast to the case of quadratic invariants, it turns out 
to be very natural to use a Kahler potential for the Jacobi metric. We find that the condition imposed 
on the geometry by the Killing tensor equations is a nonlinear PDE in the Kahler potential, whereas the 
corresponding condition for the second rank case is a linear PDE ||. Imposing the additional condition of 
invariance at arbitrary energy, the nonlinear PDE splits into a system of three equations which in general 
is overdetermined. This accounts for the fact that only isolated cases of systems with a cubic invariant at 
arbitrary energy are known. Comparing this method with the more direct approach due to Holt 0, P , we 
find the two final sets of equations to be very different. In particular, while Holt's approach leads to a PDE 
whose nonlinearity has no a priori restriction, our nonlinear PDE is always quadratic. Although our work 
therefore represents an improvement, the equations are nevertheless complicated and in practice impossible 
to solve in full generality. The fruitfulncss of our particular approach is in any case confirmed by the fact that 
we are able to present a number of new results concerning systems which are superintegrable (i.e. admitting 
the maximal number of independent invariants) at arbitrary energy. Furthermore we obtain all previously 
known strongly conserved cubic invariants, and are thus able to give a unified classification of all known 
cases. 

In addition, we have also considered strongly conserved cubic invariants of a new type which correspond 
to Killing vectors rather then Killing tensors of rank three. The existence of such "quasi-linear" invariants 
seems to be a new feature. One nontrivial superintegrable system of this type is presented, and we also 
show how its potential via coupling constant metamorphosis and coordinate translation is in fact dual to the 
harmonic oscillator potential. 
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II Jacobi geometry and Killing tensors 



The Jacobi geometrization procedure relies on the fact that given a Hamiltonian of the classical type 

H = T + V, T=\h^p aPpi V = V(q), (1) 
the orbits on a fixed energy surface H = E can be mapped onto geodesies of the Jacobi metric 

g af3 := 2(E - V)h af) (2) 



on the fixed energy surface 



Hj := y a0 Pa P = |, 9 ai 9^ = 5 a p . (3) 



The mapping is achieved via the time reparamctrization t — ► tj where dtj = 2(E — V)dt. Thus the 
Hamiltonians H and Hj represent the same system but in two different time gauges, often referred to as the 
physical and the Jacobi time gauge, respectively. For a less sketchy discussion of this time reparametrization, 
the reader is referred to Q . 

As will be made clear below a physical polynomial invariant / always corresponds to a geodesic invariant 
I j which is a homogeneous polynomial: 

Ij=K^-^ PlM --- Plim , (4) 

where K^ 1 "' IJ " m is a symmetric tensor on the configuration space, with the tensor rank m in general being 
the same as the polynomial degree of /. For the case when / is cubic, the correspondence between / and Ij 



will be given explicitly in section III , and the method used there can also be applied to polynomial invariants 
of any other degree. The equation of motion {Ij,Hj} = can readily be shown to be equivalent to the 
covariant equation 

%-^;Mm +1 ) = °> ( 5 ) 

where the parenthesis denotes symmetrization, the semi-colon denotes covariant derivation w.r.t. g a p and 
the indices of K^ 1 '"^" 1 are lowered using g a p. This is the Killing tensor equation and a solution K fll ... fim 
is known as a Killing tensor of rank m. A closely related geometric object is a conformal Killing tensor, a 
symmetric tensor -P Ml - Mm for which P fJ - 1 "' flm p f21 ■ ■ -p^ in general is invariant for null geodesies only. Since 
this is nontrivial only for the trace-free part of any tensor, a conformal Killing tensor is usually taken to be 
trace- free from the outset. The conformal Killing tensor equation then takes the form 

*W»/Wm+i) = m I 2 ( m + "] _1 .9(piM2 P,y M-M m +i);^ ( 6 ) 

where n is the dimension of the configuration space. In particular, the trace-free part of any Killing tensor 
is a conformal Killing tensor satisfying this equation. 

In this work we focus on the classical two-dimensional Hamiltonians which can be written as 

H=±(p x 2 +p y 2 ) + V(x iy ). (7) 

The associated Jacobi geometry is given by the line element 

ds 2 = 2G(dx 2 + dy 2 ) = 2Gdzdz, G = E~V, (8) 

where we have introduced the complex conjugate null variables z = x + iy, z = x — iy, which are convenient 
to use as they are adapted to the action of the conformal group. This will make the conformal Killing tensor 
equation (||) maximally simplified. Furthermore, all tensor calculations will be done in the standard null 
frame ft = G 1 / 2 dz, ft 5 = G 1 ' 2 dz in which the metric takes the simplest possible form ds 2 = 2dSl°<i£l . We 
use the convention that tensor indices in this frame take the values and 0, while in any coordinate frame 
the values will be the names of the coordinates (e.g. z and z). 



3 



Ill Mapping invariants between the physical and the Jacobi time 
gauge 

We consider a Hamiltonian of the type ([!]), admitting a cubic invariant / at least on some fixed energy 
surface H = E and possibly at arbitrary energy. The invariant can without restriction be assumed to have 
the form 

I = A a ^ PaPpPl + B a p a , (9) 

since an additional term that is even in the momenta must Poisson commute separately with H. To see how 
the invariant transforms when going to the Jacobi time gauge, we perform the time reparametrization by using 



the method of coupling constant metamorphosis 10 . To this end we wish to introduce a coupling constant 
k into the Hamiltonian by rescaling the potential according to V — ► 2kV. To see how this rescaling affects 
the invariant, it is useful to implement it in two steps, which both trivially preserve Poisson commutativity 
although they are in fact noncanonical transformations. The first step is to rescale the momenta according 
to p a —y \~ 1 p a which gives 

H -► X~ 2 T + V, 

I -> X- 3 A a ^p aP0Pl + \- 1 B a Pa , (10) 
{ , } - A{ , }, 

where { , } is the Poisson bracket. By the bilinearity of the Poisson bracket, we are also free to rescale the 
commuting functions H and / themselves. The second step is to set k = \X 2 and use this rescaling freedom 
to redefine H and I according to 

H := X 2 H = T + 2kV, (11) 
J:= A 3 / = A af3 ~<p aPl 3 Pl + 2 K B a p ai (12) 

which gives the desired rescaling of V. We could of course have defined I using any scaling factor, but A 3 
will turn out to be the natural general choice. As in Jllfl , to obtain the Jacobi Hamiltonian when solving a 
fixed energy constraint for k, we must hold on to the interpretation of the parameter E as the energy value 
of the original Hamiltonian H = T + V, thus making 2kE the corresponding energy value of the rescaled 
Hamiltonian H. Solving H — 2kE for k then results in 

K= [2(E-V)]- 1 T = H Jl (13) 

where we have used the standard definition of the Jacobi Hamiltonian Hj associated with the fixed energy 
surface of the original Hamiltonian H = T+V. Thus we find that a coupling constant metamorphosis acting 
on k is equivalent to a transformation to Jacobi time. However, as noted in pll, this is not a coupling constant 
metamorphosis in the true sense, as the old energy E does not enter linearly into the new Hamiltonian Hj. 
Nevertheless the results of j^] still apply which means that the physical invariant / and Jacobi invariant 
Ij, Poisson commuting with H and Hj respectively, are transformed into each other according to 

Ij = I u , I=Ij\ e ^ h . (14) 

Hence the Jacobi invariant corresponding to the physical invariant (||) is 

I, j = A a ^p a p pPl + 2HjB a Pa = K a ^p a p (3Pl , K a ^ := A a ^ + B^g^\ (15) 

where g af3 is the inverse of the Jacobi metric g Q/ g := 2(E — V)h a p. It follows directly that K a p^ is a third 
rank Killing tensor w.r.t. the Jacobi metric. Thus we see that a cubic invariant of the form (p) can always be 
mapped to such a Killing tensor. However, there might in fact be cases where the invariant (151) is reducible 
and for which a nontrivial cubic invariant more naturally corresponds to a Killing vector, rather then a third 
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rank Killing tensor. By appyling the transformation recipe (|l4|), one realizes that for this to happen, the 
cubic invariant (Q) must be of the form 

I = HC a Pa + D a Pa , (16) 

for some vectors C a , D a . This is clearly the same as saying that the tensor A"-^ 1 has a vanishing traceless 
part. Now, if after taking p a — > X~ 1 p a we would have defined the rescaled invariant ( |l6| ) according to I := XI 
instead of I — * A 3 / as was done in eq. (|l2|), the corresponding Jacobi invariant would simply become 

Ij = CPa, C=EC a +D a (17) 



from which it follows directly that £ a is a Killing vector w.r.t. the Jacobi metric. In section VI we give an 
example of a Hamiltonian which has a nontrivial cubic invariant of the type ( |l6| ) . 

IV Cubic invariants corresponding to third rank Killing tensors 

In this section we begin by deriving a necessary and sufficient integrability condition for the Jacobi metric to 
admit a third rank Killing tensor at a fixed value of the energy parameter E. In the next section we proceed 
by finding the conditions that ensure that the Killing tensor equations will be satisfied at arbitrary energy 
values, but with attention restricted to the case when the energy dependence of the Killing tensor is such 
that the corresponding invariant of the physical Hamiltonian is cubic. In analogy with the second rank case 
H as well as the third rank case with indefinite metric (l2) , our Killing tensor K a p 7 will from the outset be 
decomposed into its trace-free (conformal) part P a p~ t and trace K a := K@ p a according to 

K af 3 7 = P Q/ 3 7 + 3 (n + 2y 1 K( a g l3l ), (18) 

where n is the dimension of the configuration space, i.e. n — 2 in our case. Similarly the Killing tensor 
equation itself, Kt a p Ti $) — 0, will be split into its three trace-free "components", namely its trace-free part, 
the trace-free part of its trace and the trace of its trace, which read 

C Q( 3 7 5 := P( a p T ,8) - 3 (n + 4) _1 .g (ct/3 F A 75) . A = (19) 
D aP := K {a . 0) + (n + 2) [n + 4t)~ 1 P r ' a p ri = (20) 
K a , a = 0, (21) 



where the trace of eq. (|20|) is automatically satisfied once the divergence-free condition (|21j) for K a has 
been solved. Just as in the case of an indefinite metric |1^ |, it is advantageous to use the coordinate 
frame components of K a and P q ' 3 t (note the index positioning) when parametrizing the four independent 
components of K a p 7 . The parametrization thus becomes 



Kooo = G 3 / 2 P^ Z 
K m - G 3 / 2 P^ 

V V-l/2 TS 

iv 000 — 2 U z 



(22) 



The calculations leading to the integrability condition imposed on the Jacobi metric are analogous to the 
indefinite case |l2| ]. An important difference, however, is the fact that the null variables z and z here are 
complex conjugate. Demanding that the potential V (and thereby the Jacobi metric) be real implies that 
there is no restriction in assuming that the Killing tensor K a py be real as well. This in turn leads to the 
component constraints that Kg and P zzz be the complex conjugates of K z and P ZZZ J respectively. In the 
indefinite case, on the other hand, the reality condition for the Killing tensor requires that the corresponding 
four components be real, but implies no other relation between them. 



5 



We now proceed by solving the two equations ( |l9| ) and (|T]), leaving (^Oj) - the only equation which 
couples the conformal and trace parts of the Killing tensor - for later consideration. 
The conformal Killing tensor equation ( |l9| ) for P Q /3 7 has the components 



Coooo — Coooo — GP ZZZ ^ — 0. (23) 

Hence P a /3y can be fully represented in terms of an analytic function S(z) :— p zzz . We shall assume S(z) 7^ 
to be the case, since otherwise eq. (|o|) implies that K a is a Killing vector. The divergence free condition 
(0 for K a , 

K Q . a = 2G- 1 Re{K z ^} = Q, (24) 

is just as easily solved by introducing a real potential function $ satisfying K z — 2i$ iZ , with the factor 
2 inserted for later convenience. This potential relation can also be expressed covariantly in terms of the 
natural volume 2-form e a p = A f2°) Q/ 3 as 

K a = 2e/$ ;/3 , (25) 
provided that <!> transforms as a scalar. We can now write eq. ( p0| ) in the form 

D al3 = 2$. 7(Qe/3) T + §P 7 Q/3:7 = 0. (26) 

Its components, 



Dun = Z) (11) = -2i(^), z - + |G- 2 (G 3 5),, = 0, (27) 



can in analogy with the procedure used in J12[ be simplified by using a conformal transformation to a new 
complex null variable w — H(z), with the inverse relation being z = F(w). The metric will then be given 
in terms of the new conformal factor G = \F'(w)\ 2 G as ds 2 = 2Gdwdw. To preserve the relation between 
the standard null frame and the null variable, the frame must be scaled as Ct° = Bfl°, f2° = B^ 1 ^ where 
B^ 1 = B = [F' '(w) J 'F '(w)] 1 / 2 . The analytic function S(z) transforms as 

S(u) := P www = [H'{z)fP zzz = [H'(z)] 3 S(z). (28) 

This shows that we can always make S(w) take the standard constant value 1 by choosing the conformal 
transformation such that H'(z) = [S(z)]^ 1 ^ 3 . With this choice, eq. (27) simplifies to 

A)0 = -2i(^),»+2G, w = 0. (29) 
Comparing this equation with its complex conjugate leads directly to the integrability condition 

M(%),—} = 0, (30) 

which is satisfied iff $ and G are related by some real function JC according to 

$ U! = iGK^w. (31) 

Substituting this back into eq. ( |29| ) yields 

G _w — fc,www- (32) 

In fact, as K, is determined only up to the transformation JC — * K, + Avow + Re{A(w)}, where A is a real 
constant and A(w) is an analytic function, we can partially fix this freedom by requiring that the relation 

G = JC^ W w (33) 



G 



holds, meaning that K. becomes a Kahler potential |T3|] for the Jacobi metric. As such, JC transforms as 
a scalar under conformal transformations. Hence for any null variable z, the metric conformal factor G is 
related to K, according to the simple formula /C z5 = G. Substituting eq. ( |33| ) into eq. ( |3l| ) now leads to the 
final integrability condition 

Re{(lC iWiD }C lWW ) tW } = 0, (34) 

which is necessary and sufficient for the existence of a third rank Killing tensor. In contrast to the first and 
second rank cases ||, we are here dealing with a nonlinear condition imposed on the metric, ft can be noted 
that a condition which is formally identical to eq. (|^) also applies to one of the three different types of third 
rank Killing tensors that are allowed when the metric has indefinite signature |Q. The interested reader 
can consult that reference for a number of nontrivial solutions. Moreover, the same condition can be found 
in H , but not in the context of the geometric approach of this work in which the unknown function can be 
interpreted as a Kahler potential. 

Going back to the original null variable z, eq. ( |3l| ) now takes the form 

= t/C, rf (S/C,« + ±S7C, 2 ), (35) 

with the corresponding integrability condition 

Re{[AC,«(SAC,« + |S"fc,,)],,} = 0. (36) 

The cubic Jacobi invariant Ij — K a ^ 1 p a ppp 1 can now, via eq. ([35]) and G = lC tZZ , be written down as 

Ij = 2 Re{S Pz 3 - {3S1C tZZ + S'K z )H jPz }, Hj = K^PzPz- (37) 

Expressed in the standardized null variable w, this simplifies to 

Ij = 2Re{p„, 3 - 3K tWW Hjp w }, Hj = IC^w^PwPw- (38) 

Note that up to conformal transformations, the Jacobi Hamiltonian and its cubic invariant is completely 
determined by K. 



V Arbitrary energy invariants 

So far we have derived the necessary and sufficient condition (eq. ( |34"| ) or eq. (|36|)) for a two-dimensional 
Riemannian geometry to admit a third rank Killing tensor. Given that this condition is satisfied, we can use 
any null variable z and identify the associated metric conformal factor G — ]C jZZ with E — V for a flat space 
Hamiltonian H — 2p z p z + V which is integrable with cubic invariant at fixed energy E. Thus we can choose 
to interpret the geometry as a common Jacobi geometry of a large family of such flat space Hamiltonians, 
whose members are related by conformal transformations. Now, however, we shall make the identification 
^C,z2 = E — V for a particular null variable z and derive the additional conditions that make the Killing 
tensor equations satisfied for all values of E, corresponding to H = 2p z p s + V being integrable at arbitrary 
energy. To this end we begin by noting that we can write the Kahler potential as 

JC = Ezz + 2Rc{9(E,z)} - (39) 

where 6(E, z) is analytic in z and can be assumed to satisfy 9(0, z) — 0, whereas VP is real, satisfies ^ :ZZ - = V 
and is independent of E. It appears at first sight that 0(E, z) in principle could have any dependence on 
E, which would make it very difficult to continue working in full generality. Similarly it seems possible that 
the function S(z) could also have some dependence on E. However, if we apply the recipe that the Jacobi 
invariant Ij is transformed to the physical time gauge according to I = Ij\e^h^ ^ follows from eq. ( |36| ) 
and eq. (^) that / becomes a cubic polynomial invariant iff S(z) is independent of E and no restriction is 
implied by setting 6(E, z) — EA(z) for some analytic function A(z) which is also independent of E. Since 
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this paper is devoted to cubic invariants, we shall use this fact from the outset. Our physical cubic invariant 
/ then has the general form 

I = Re{2S Pz 3 - 2[S'(z + A') + 3SA"] Pz 2 P - z + - (z + A')9 ltg ) + 3S(*, ZZ - A"* >zz )]p z }. (40) 

If we now substitute 

JC = E[zz + 2Re{A{z)}] - * (41) 

into eq. (|3^), the condition becomes a second degree polynomial in E: 

A 2 E 2 + A X E + A = 0, (42) 

where the coefficients A^ must vanish separately if the equation is to hold for arbitrary E. This gives the 
three equations 

A 2 = Re{[5'(z + A') + 35A"],4 = 0, (43) 

A 1 = - Re {[£'((« + A')*,,* + ¥,,) + 3S(A"¥, a * + *, zz )], z } = 0, (44) 

A Q = Re{[^>MS'^,z + 3Stf ,„)],*} = 0. (45) 

By applying the differential operator d 2 jdzdz to eq. ([43]), we split off the condition 

Re{S"'(z)} = 0, (46) 

with the polynomial solution 

S(z) =iaz 3 + f3z 2 +jz + S, (47) 

where a is a real constant and (3, 7 and S are complex constants. A standardization of these coefficients can 
be achieved by using the transformation property of S(z) given by eq. (p8|), combined with the freedom to 
make translations z — > z + Zq, rotations z ~ * e lb z and scalings z — > cz, with & and c real. Comparing with 
the first and second rank cases B leads us to conjecture that the analytic function S(z) representing the 
conformal part of a Killing tensor of arbitrary rank m is required to satisfy 

Re{(±rS(z)} = (modd), 
Tm{(£) m S(z)} = (meven). 

However, we make no attempt to prove this. For any choice of coefficients a, (3, 7 and 5, eq. (^) is a linear 
ODE in A(z) which can be solved by standard methods. Once this first equation has been taken care of, eq. 
( fli] ) and (45) are two PDEs for one unknown function ^P, i.e. in general an overdetermined system. This 
is consistent with the fact that only isolated cases of classical Hamiltonians integrable with cubic invariant 
are known B . The obvious way to proceed is to try to find the general solution to the linear PDE ( |44| ) and 
then check if its functional degrees of freedom can be exploited to make the nonlinear PDE ( Eg) satisfied. 
As a simple illustration of this method, we start out from the simplest possible solution to eq. (|43|), namely 
a = /3 = 7 = A(z) = 0, 8 = d with d real. In terms of the real variables x and y which satisfy z = x + iy, the 
linear condition ([l4j) in this case becomes 

" ^,xyy = 0. (49) 

Writing down the general solution as 

* = h(y + V3x) + f 2 {y - V3x) + f 3 (-2y), (50) 
the nonlinear condition ( f45"l ) takes the cyclic form 

F X '(F 2 - F 3 ) + F 2 '{F 3 - Fi) + F 3 '(Fi - F 2 ) = 0, (51) 
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where Fi — //' so that 

V = F 1 {y + V3x) + F 2 {y - V3x) + F 3 (-2y). (52) 

The functional condition is well-known from Lax pair studies ^ |l4| and arises for a three-particle 
Hamiltonian of the generic form 

H = \W + P2 2 +P3 2 ) + Fxiq 3 - q 2 ) + F^q 1 - q 3 ) + F 3 (q 2 - q 1 ), (53) 



which can be reduced to our two-dimensional Hamiltonian with V given by eq. ( p2| ) by going to the center- 
of-mass system and making a rescaling of the coordinates ||. Its solutions include the Toda potential 
Fi(£) = Cie H and also = CV(b£), where the Weierstrass function V(£) can be taken in the special 

limits — > £~ 2 , T^iO — * sinh" 2 £ and V(^) — > sin~ 2 £. The cubic invariant of a system of this type can 
be found in table [l]. 

In more general cases, the linear PDE (|44| ) is not as easily solved. Failing to work in full generality, we 
have restricted attention to the special cases where S(z) is a homogeneous polynomial, i.e. we have set all but 
one of the constants a, /3, 7 and S to zero. Furthermore, for all our solutions, whether previously known or 
not, it has turned out that A(z) = Xz 2 for some complex constant A. Thus the coefficient zz + 2 Re{A(z)} of 
E in the Kahler potential ( fid]) is always a hermitian form in z (or equivalently, a real quadratic form in x and 
y). However, we stress that solutions of other types may very well exist. With both S(z) and zz + 2 Re{A(z)} 
being homogeneous functions of z, z, it is very natural to make the anzats that \& is homogeneous as well. 
To this end one can e.g. introduce polar coordinates r, <f> according to z = re 1 ^ and set 'J = r k /(</>), where 
k is a real constant and f{4>) is an arbitrary real function. This form of VP gives a Jacobi metric 

ds 2 = 2(E - V){dr 2 + r 2 d(f> 2 ), V = \r k - 2 [k 2 f(</>) + /"(<£)] (54) 

which, at zero energy (all energies iff k = 2), admits a homothctic vector field C = 2fc~ 1 r<9/9r, satisfying 
£(9a/3 = 2g a {3. The advantage of this ansatz is that eq. ( f44| ) and ([f5j) become ODEs in f(<fi), parametrized 
by k. Several of our solutions, some of which are new, were found using this approach. However, it is worth 
mentioning that in some cases the ODEs become simpler in other coordinates (such as £ = z, r\ — z/z) which 
are also adapted to the homogeneous ansatz. 

VI Cubic invariants corresponding to Killing vectors 

In this section we look for geom etrie s admitting a Killing vector £ a corresponding to a cubic invariant by the 



mechanism discussed in section III. If we parametrize the Killing vector by £ 2 and its complex conjugate, 
the trace-free part of the Killing vector equations = (i.e. the conformal Killing vector equations) 

becomes 



£(0;0) = £(0;0) = C",S = 0, (55) 

which is solved by setting £ z = Z(z) with Z(z) analytic. The linear Jacobi invariant thus takes the form 
I J = i a Pa = 2Kc{Z(z)p z }. The trace part of the Killing vector equations becomes 

= 2G- 1 Re{(S*G) :Z } = 0. (56) 

Under conformal transformations, S(z) transforms according to 

Z(w) := C = H'(z)^ z = H'(z)Z(z). (57) 

Hence we obtain Z(w) = 1 by choosing the transformation so that H'(z) = [Z(z)]^ 1 , which makes eq. (^) 
take the standardized form 

Re{G = 0. (58) 
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It follows that G is a function of Y = Ivaw only, which is consistent with the invariant being Ij — 2px, 
X = Rew. 

We now fix a null variable z for which we set G = E — V and focus on trying to make the Killing vector 
equations hold for arbitrary energy E, for the case that the Killing vector corresponds to a cubic physical 
invariant. According to eq. ( )l7j ) we must allow for Z(z) to have a linear dependence on E: 

Z(z)^ Z(E,z) = Z 1 (z)E + Z (z), (59) 

where we have set Z\{z) = C z and Zq(z) = D z . This leads to the physical invariant taking the general form 

I = 2 Re{(Z x H + Z )p z }, H = 2 PzPz + V. (60) 

Substituting eq. ( |59] ) into eq. ( |56| ) yields 

iRe{(ZG), z } = B 2 E 2 + J B 1 £ + Bo = 0, (61) 



giving us the three equations 



The solution to eq. (p2h reads 



B 2 = RelZ/} = 0, (62) 
B 1 = Re{(Z - Z t V), z } = 0, (63) 
B = -Re{(Z V), z } = Q. (64) 



Z\ [z) = iaz + j3, (65) 



where a is a real constant and (3 is a complex constant. This gives us two inequivalent cases to study since 
we can make a linear transformation such that cither Z\{£) — * 1 or Z\{z) — > iz depending on whether a is 
zero or nonzero. We consider these two cases in some detail below. 

The Case Z t (z) = 1 

In this case eq. (|63| ) takes the form 

Rc{Z ' - V z } = 0, (66) 

with the general solution given by 

V = Z (z) + Z (z) + F(y), y = lmz, (67) 
Here F(y) is an arbitrary real function. The final condition ( [64] ) now becomes the functional equation 

Re{[Z (Z + Z a + F)], z } = {Z + \Z + \F)Z Q ' + (Z Q + \Z + \F)Z ' - \(Z Q - Z )F' = 0. (68) 
We have found no nontrivial solution to this equation. 

The Case Z x {z) = iz 



Here eq. (44) reads 

Re{Z ' - izV z } = 0, (69) 



and is solved by 



V = i I - / — dz + f ^-dz ) + F(r), r = ^zE, (70) 



z z 
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with F(r) being an arbitrary real function. The integrals can be avoided by introducing an analytic function 
Q(z) satisfying 



Z = i(zQ'-Q). (71) 

We can then write the potential as 

V = Q'(z) + Q'(z) + F(r), (72) 

which makes eq. (p4) become 



Im{[(zQ' - Q)(Q' + Q' + F)], z } = ±(2zQ' - Q + zQ' + zF)Q" - ^(2zQ' - Q + zQ' + zF)Q l 

+ U r Q' ~ r Q' ~ zQ/ r + *Q/r)F' = o. 



(73) 



This functional equation has at least one nontrivial solution, given by Q(z) = a^/z, F(r) = Ar 1 with a 
complex and A real. Writing a = V2(B + iC), the potential becomes 



v = A + BVFT^ + CVT^lc 
r 

This potential was given in complex form by Drach ||. We also note that the potential also admits a 
quadratic invariant |5j and is therefore now superintegrable. The cubic invariant of the system, given by 
cq. (|6(]) , is presented in table |[ What is interesting about this result is that it gives a new way of solving 
Hamilton's equations for the system, even in the physical time gauge. Namely, for any fixed energy E, 
one can choose coordinates such that the cubic invariant / reduces to the momenta pq of a configuration 
coordinate Q that becomes cyclic in the Hamiltonian for this particular energy value. In other words the 
equations of motion can be solved just as in the case of an ordinary linear invariant, but the coordinate 
transformation which makes Q cyclic is dependent on the energy E. In the geometric picture, one obtains 
a cyclic variable by the above described conformal transformation which makes the metric conformal factor 
satisfy eq. (p^) . 

Quite surprisingly, we soon realized that this potential via coupling constant metamorphosis acting on 
A is dual to the usual (isotropic) harmonic oscillator potential, with its center in general translated off the 
coordinate origin. This coupling constant metamorphosis can be realized as a conformal transformation in 
the Jacobi geometry setting. Choosing the transformation as z — F(w) — w 2 , the metric conformal factor 
G = E — V with V given by eq. ((TJ) , transforms into 

G = \F'\ 2 G = E- V, V = A(X 2 +Y 2 ) + BX + CY. (75) 

where X = Re w, Y = Im w, E = -4A, A = -4E, B = 4y/2B and C = 4V^C. We have assumed that X and 
Y are both positive to avoid keeping track of signs. It follows that V can be interpreted as a potential which 
is integrable at arbitrary energy E. Obviously this V is just the translated harmonic oscillator potential. 

The analytic function Z(E,z) = i[Ez — 2~ 1 / 2 (B + iC)^/z\, on the other hand, transforms according to 
eq. ( |57| ) into 

Z{w) = —^(2Aw + B + iC). (76) 
Since Z(w) is independent of the new energy E, the invariant is the same in both time gauges and reads 

I = Ij= M 2Ay + C)P* - ( 2Ax + B)Py], (77) 

where, of course, the factor jg could be dropped by rescaling Z(w). Clearly / is the ordinary, well-known 
linear invariant for the system B. 
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VII Classification of Hamiltonians admitting a cubic invariant at 
arbitrary energy 

In this section we present our nontrivial solutions to the equations ( fi"3| ) - ( f45| ) and (^) - ( |64[ ) . For the third 
rank Killing tensor cases the solutions can be fully represented by the analytic function S(z) and the Kahler 
potential /C. The solutions are given in four tables (|-|), one for each possible degree of the polynomial S(z). 
For the single Killing vector case, the analytic function Z(E, z) and the potential V give all information. 
This solution is given in table M. 

For all cases we will present the potential V (up to linear transformations of the coordinates and a 
rescaling of the potential itself) as well as the physical cubic invariant /, thereby making it possible to 
directly compare our results with Hietarinta's classification of 1982 ||. A real scaling of S(z) or Z(E,z) 
only results in an irrelevant scaling of the invariant / by the same factor. Accordingly this freedom will be 
fixed such that I takes a convenient form. 

For the cases which are superintegrable with both quadratic and cubic invariant we also indicate, using 
the notation of which of the real quadratic cases (1), (2), (4) or (7) a given system belongs to. 

VIII Comments 

We have shown that cubic invariants at fixed and arbitrary energy can be treated in a unified manner by 
using the Jacobi geometrization method. Most strongly conserved cubic invariants are nontrivially cubic in 
our geometric picture, but we also found a mechanism by which a new type of cubic invariant can instead 
correspond to a family of linear invariants parametrized by the energy. It is then possible to obtain a cyclic 
variable by making a standardizing conformal transformation which is energy dependent. A natural first 
extension of this result could be to investigate if the same mechanism can be explored to find nontrivial 
quartic invariants which analogously correspond to quadratic invariants standardized by means of energy 
dependent conformal transformations. This is particularly interesting considering that the standardizing 
conformal transformations in the quadratic case are associated with explicit Hamiltonian separability 
Compare also recent work by Rauch-Wojciechowski and Tsiganov [jl7j who gave some examples of non- 
standard separability. In fact they considered a more general situation with the separating transformation 
involving the second invariant itself in addition to the energy. 

It would also be of interest to apply the method used in this paper to the usual type of quartic invariants 
which does not reduce to quadratic invariants when fixing the energy. In [[l8) it was shown that the integra- 
bility condition for fourth rank Killing tensors is of the same nonlinear type as in the third rank case. Thus 
it should be possible to impose the arbitrary energy condition using the approach adopted in the present 
work. 
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S(z) = 4 




JC = E(x 2 + v 2 ) - h(v + VSx) - f 2 (y - VSx) - f 3 (-2v), Fi(£) := f/'(0 G |e« 


V(£), £~ 2 , sinh _2 £, sin _2 fi 


V = Fi f w + a/3 art + F9 f v - a/3 + F-i(-2v) 




/ = Pr 3 - 3»^»„ 2 + 3LFi fw + %/3aO + Fo(v - \/3x) - 2F-i(-2v)h^ + 3a/3[-Fi fw + 


a/3 x} -(- Fo( v — a/3 xWdi, 


= £ -2 is superintegrable, case (2) 




= 1 + i 




JC = E(x - y) 2 - i|(x 5 / 2 ± v 5 / 2 ) 




V = \/x± */v 




I = P-r 3 — Vi? + 3(\[xV x T ^/VVn) 
f x f y i " V V r * l \ v r y 1 




Superintegrable, case (7) 




S(z) = 1 




JC = 2Ey 2 - i|x 5 / 2 - 28~ 1 xy - §<5y 3 

y 15 y 3 y 




V = \fx + Sv 




J = rj 3 + 3\/ot^ — ^Di, 




Superintegrable, case (7) 




S(z) = -1 




AC = |F(2a; 2 + y 2 ) - i(8x 4 + 2Ax 2 y 2 - y 4 ) + 451n|y| 




V = Ax 2 + y 2 + Sy- 2 (Drach 1935 §, Holt 1982 §) 




I = PxPy 2 + 2(-y 2 + 5y~ 2 )p x + 8xyp y 




Superintegrable, case (4) and (7) 




S(z) = -1 




JC = \E{2x 2 + y 2 ) - |(|a; 3 + xy 2 ) + A5\n\y\ 




V = x + Sy 2 (Drach 1935 [[|, can be linearly combined with the potential above) 


I = PxPy 2 + 2Sy 2 p x + yp y 




Superintegrable, case (4) and (7) 




S(z) = -1 




JC = 2E(2x 2 - y 2 ) + f| y 10 / 3 - 9{x 2 + % 4 / 3 




V = |y 4 / 3 + (x 2 + 6)y- 2 / 3 (Drach 1935 @, Holt 1982 |]) 




/ = 2p, 3 + 3p^ 2 + 3[-3y 4 ^ + 2(x 2 + 5)y- 2 ' 3 ]p x + l%xy l / 3 p y 




S(z) = -1 




JC = 2E(2x 2 - y 2 ) - 9xy~^ 3 




V = xy~ 2 l 3 (Drach 1935 ||, can be linearly combined with the potential above) 




I = 2p x 3 + 3p xP y 2 + Qxy- 2 ' 3 p x + 9y^ 3 p y 





Table 1: Systems for which S(z) is of the zeroth degree. 
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S(z) = —iz 

K = \E(3x 2 + y 2 ) - §x 4 - \x 2 y 2 - ^ 

V = 9x 2 + y 2 (Fokas and Lagerstrom 1980 
I = (xp y - yp x )py 2 + lv z p x - Qxy 2 p y 
Superintcgrable, case (7) 

S(z) = -z 

K = E(x 2 + y 2 ) ~ ±A(x 4 + Ax 2 y 2 + ?/) + 4(5 ln\x\ + Cln|y|) 

V = A(x 2 +y 2 ) + Bx- 2 + Cy- 2 (Drach 1935 §) 

J = [xpy - yp x )(p x p v + 2Axy) - 2Byx~ 2 p y + 2Cxy~ 2 p x 
Superintegrable, case (2) and (7) 
S(z) = 2iz 

K = Er 2 - 6 |sin(2</>/3) / cos(20/3)[cos(20)]- 2 / 3 d</> - cos(2<£/3) / sin(20/3)[cos(20)]- 2 / 3 d</>} r 2 / 3 

V = (x 2 - y 2 )~ 2/3 (Fokas and Lagerstrom 1980 @) 

I = (xp y - yp x ){p x 2 - Py 2 ) - A(yp x + xp y ){x 2 - y 2 )~ 2/3 

S(z) = 2iz 

K = E(x 2 + y 2 ) - ±A(3x 4 + 2x 2 y 2 + 3?/) - ±Bln| (x + y) (x - y)- 1 ] + Cm \x 2 - y 2 \ 

V = A(x 2 + y 2 ) + [Bxy + C(x 2 + y 2 )](x 2 - y 2 y 2 

I = (xp v - yp x )(p x 2 - p y 2 ) + 2A(xp y - yp x )(x 2 - y 2 ) 

-[x(x 2 + 3y 2 )(B Px + 2Cp y ) + y(3x 2 + y 2 )(2C Px + Bp y )]{x 2 - y 2 y 2 
Superintegrable, case (2) 



Table 2: Systems for which S(z) is of the first degree. 



S(z) = -z 2 




JC = iE(2x 2 + 3y 2 ) - A[Ar - B ln(r + x) - C ln(r - x)\ 




V = [A + B(r + x)- 1 + C(r - x)" 1 ]^ 1 




I = (x Py - yp x ) 2 Px + 2{B + C)[l + (x/y) 2 ] Px + {-Ay(xp v - yp x ) + (C - B) ( [2(x/y) 2 + 3]xp x - 


Vyp y )}r 1 


Superintegrable, case (2) and (4) 




S(z) = -z 2 




K = E[l - \ cos(2</>)]r 2 - i(Jf((/)) dcf)-C lnr), 




condition: [3/"/' - 2/'/ - Cf"} sin + [/"/ + 4(/') 2 - 2Cf] cos = 




V = f^)r- 2 




I = p 2 (cos(/)p r - sin^r -1 ^) + [2/'(» cos - f(<f>) smcj)]p r - [3/'(0) sin0 + f(4>) cos^r -1 ^ 




Superintegrable, case (2) 




S{z) = -z 2 




K = \E{2x 2 + 3y 2 ) - A(Ar + B In \ y (r — — C In \y ) (special case of the above when A 


0) 


V = Ar- 1 + (Bxr- 1 + C)y- 2 (Drach 1935 §]) 




/ = (xp y - yp x ) 2 p x - A(xp y - yp x )yr- x + [B(2x 2 + 3y 2 )xr" 1 + 2Cr 2 ]y- 2 p x + Byr^py 




Superintegrable, case (2) 





Table 3: Systems for which S(z) is of the second degree. 



14 



S(z) = -iz 3 

K, = Er 2 — Af{(j>)r-\ condition: /"'/" - 2f"f - 3/'/ = 

V = [f{4>) + f"((t))}r- 3 (compare with Thompson 1984 JT|) 

/=^ 3 + 3TO)p r + r(0)r-y ] 

S{z) = -iz 3 

JC = Er 2 - A[{\nr) 2 + cj) 2 ] - (Be^ + Ce'^y- 1 (special case of the above when A = 0) 

V = Ar- 2 + (Be^ + Ce-^y- 3 (Drach 1935 §, Hietarinta 1986 [§) 
I = P0 3 + 2^(Be^ - Ce-^) Pr + \[2A + 3(Be^ + Ce-^y- 1 ^ 

Table 4: Systems for which S(z) is of the third degree. 

Z(E, z) = i[Ez - 2- 1 / 2 {B + iC)sfz\ 

V = {A + B^/TT^+C^/T^y- 1 (Drach 1935 @) 

/ = (xpy - y Px )H + \ [(By/r^x + Cy/7T^)p x + {-B^/r + x + Cy/7^x~)p y ] 
Superintegrable, case (4) 

Table 5: A system with a cubic invariant corresponding to a Killing vector. 
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